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Abstract 



Q ' We consider random hermitian matrices made of complex blocks. 

CJ ■ The symmetries of these matrices force them to have pairs of oppo- 

site real eigenvalues, so that the average density of eigenvalues must 
vanish at the origin. These densities are studied for finite N x N 
j_j ■ matrices in the Gaussian ensemble. In the large limit the density 



of eigenvalues is given by a semi-circle law. However, near the origin 
there is a region of size in which this density rises from zero to the 
semi-circle, going through an oscillatory behavior. This cross-over is 
calculated explicitly by various techniques. We then show to first or- 
der in the non-Gaussian character of the probability distribution that 
this oscillatory behavior is universal, i.e. independent of the proba- 
bility distribution. We conjecture that this universality holds to all 
orders. We then extend our consideration to the more complicated 
block matrices which arise from lattices of matrices considered in our 
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previous work. Finally, we study the case of random real symmetric 
matrices made of blocks. By using a remarkable identity we are able 
to determine the oscillatory behavior in this case also. The universal 

oscillations studied here may be applicable to the problem of a particle 
propagating on a lattice with random magnetic flux. 
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1 Introduction 



It is well known that the average density of states, for Gaussian ensembles 
of random N x N matrices, obeys Wigner's semi-circle law when goes to 
infinity, irrespective of the symmetries of the probability measure For 
no n- Gaussian measures the average density depends sensitively upon the 
distribution [Q. However, next to the edge of the support of the eigenvalue 
distribution, there is a region of size N^'^^^, in which the average density 
crosses over from nonzero to zero, with a universal cross-over function (i.e. 
independent of the probability distribution) [Q. In this work we consider 
instead an ensemble of random hermitian matrices made of complex blocks. 
These matrices have been discussed recently for its application to impurity 
scattering in the presence of a magnetic field 0, |^ and to a study of the 
zero modes of a Dirac operator 0. In the large N limit the average density 
of eigenvalues is again a semi-circle for Gaussian ensembles. However, by 
construction these matrices have pairs of opposite real eigenvalues. Thus, as 
an eigenvalue approaches zero, the mid-point of the spectrum, it is repelled 
by its mirror image. Consequently, the density of eigenvalues is constrained 
to vanish at the origin. Away from the origin, the density must rise rapidly, 
over a region of size towards the Wigner semi-circle. 

In a recent work we showed, by explicit computation in the Gaussian 
case, that the rise "overshoots" the Wigner semi-circle and thus has to come 
back down, whereupon it overshoots again. Thus, the density of eigenvalues 
oscillates over a region of size . This cross-over at the center of the spectrum 
is however not of the same nature as the cross-over at the end of the spectrum. 

For the simplest case of one random matrix, we can calculate the cross- 
over by three different methods: i) the orthogonal polynomial approach [0], 
ii) a method inspired by Kazakov's approach to the usual hermitian Gaussian 
problem , iii) a supersymmetric method based on Grassmannian variables 
[^]. The first method is in fact quite cumbersome and hard to generalize to 
more difficult problems. The second one provides an elegant integral repre- 
sentation of the correlation functions for finite A^. However it is only through 
the Grassmannian approach that we could tackle the more complicated prob- 
lems of a whole lattice of matrices discussed in our previous work ||19|| . 

For the most part the calculations in this paper are done with the Gaus- 
sian distribution. An interesting question is whether the cross-over behavior 
at the center of the spectrum is universal, that is, independent of the details 
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of the probability distribution. A simple scaling argument suggests univer- 
sality. As usual in random matrix theory , the eigenvalues can be thought 
of as representing gas particles in a one-dimensional space. The probability 
distribution of the random matrices determine the potential confining the 
gas. This confining potential controls the width of the eigenvalue spectrum, 
which is of order in our convention. Crudely speaking, the oscillatory be- 
havior we are discussing here depends on the repulsion between an eigenvalue 
and its mirror image, and on the repulsion between the eigenvalues. It should 
be possible to regard the confining potential as essentially constant over the 
region of size -j^ we are concerned with and hence irrelevant. However, fur- 
ther thought reveals that this argument is insufficient, since the confining 
potential, by changing the width of the spectrum, effectively also changes 
the average value of the density of eigenvalues near the center of the spec- 
trum. We also expect that as the width of the entire spectrum changes the 
period of the oscillations near the center of the spectrum would also change 
accordingly. We are able to show, to first order in the non-Gaussian charac- 
ter of the distribution, that these effects cancel out, and that the cross-over 
behavior is universal. We conjecture that this cross-over behavior is indeed 
universal to all orders. 

We then show that we can extend our considerations to study random 
real symmetric matrices made of real (but not symmetric) blocks. By using a 
remarkable identity, the universal oscillations near the center of the spectrum 
can again be calculated explicitly. 

Finally, we discuss a possible realization of these universal oscillations in 
the problem of a single particle propagating on a square lattice penetrated 
by a random magnetic flux. 

2 Orthogonal polynomials approach 

In the simple one-matrix problem we consider 2N x 27V block matrices M of 
the form 



in which C is an x complex random matrix, with probability distribution 




(2.1) 




(2.2) 
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It is easy to show that M has pairs of opposite real eigenvalues; indeed if 
^1 is an eigenvector for the eigenvalue A, then ( ^ ) is an eigenvector 

yJ \-yJ 

for —A. In other words the matrix M anti- commutes with the "75" matrix 
^ ^ . Therefore one can express the average resolvent of the matrix M 



-1 

in terms of that of C^C, a hermitian matrix with positive eigenvalues: 



Taking the imaginary part of ( p.3|) we relate the density of eigenvalues of M 



to that of etc 



P(A) =< ^Tr5(A - M) > (2.4) 



p(r) =< ^Tr5(r - C^C) > (2.5) 



as 



p(A) = \X\p{X'). (2.6) 

If we integrate out as usual over the unitary group, we are left with 
integrals over the N eigenvalues of C^C, which run from zero to infinity. 
The Jacobian of this change of variables is simply the square of the Van der 
Monde determinant of the r^'s. The probability measure is then 

1 ^ 

P{ri, ■ ■ ■ , r;v) = -exp(-iV5]ri) A2(ri, ■ ■ ■ , r^). (2.7) 

^ 1 

All the correlation functions are known to be expressible in terms of the 
kernel 

K{r, s) = exp( ^^) ^ L4Nr)L^iNs). (2.8) 



in which the polynomials Ln{x) are the Laguerre polynomials, orthogonal on 
the half-line with the measure exp(— x). The density of eigenvalues of C^C 
is then given by 

< ^Tr(5(r - C^C) >= K{r, r) . (2.9) 



5 



Using the Christoffel-Darboux identity, and the asymptotic behavior of 
Ln{Nr) for n of order N, we could obtain the desired density of state from 
this method. But not only is this method more cumbersome for this one- 
matrix case, also it cannot be generalized to the lattice of matrices that we 
discuss below. We shall return to it later after we considered the cross-over 
near the origin. 

3 Kazakov's method extended to complex ma- 
trices 

3.1 Contour integral 

For the usual Gaussian ensemble of random hermitian matrices, Kazakov 
has introduced a curious, but very poweful, method ||^ . It consists of adding 
to the probability distribution a matrix source, which will be set to zero at 
the end of the calculation, leaving us with a simple integral representation 
for finite A^. However one cannot let the source go to zero before one reaches 
the final step. Let us follow the same lines for the model at hand. We modify 
the probability distribution of the matrix by a source A, an N x N hermitian 
matrix with eigenvalues (ai, • • ■ , cat) : 

P^(C) = -^expi-NTiC^C - NTtAC^C). (3.1) 

Next we introduce the Fourier transform of the average resolvent with this 
modified distribution: 

UAit) =< ^Tre^*^'^ > (3.2) 
from which we recover, after letting the source A go to zero, 

/OO fif 
_^^^"''Uo{t). (3.3) 

Without loss of generality we can assume that A is a diagonal matrix. We 
first integrate over the unitary matrix V which diagonalizes C^C . This is 
done through the well-known Itzykson-Zuber integral over the unitary group 
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/--(T'-^"')-^^ (3.4, 

where A (A) is the Van der Monde determinant constructed with the eigen- 
values of A: 

A{A) = l[{a,-aj). (3.5) 

i<j 

We are then led to 

1 1 ^ r 

= z:a(I)]v£/*'---*"^*"^''-'---''"> 

N N 

X exp{-NY.r,-Nj2rjaj). (3.6) 

i=i j=i 

Then we have to integrate over the rj's. It is easy to prove that 

N 

/ dri--- drNA{ri, ■ ■ ■ , tat) exp(- ^ r^^i) 

1=1 

With the normalization Ua{0) = 1, we could always divide, at any interme- 
diate step of the calculation, the expression we obtain for UA{t) by its value 
at t = 0, and thus the overall multiplicative factors in (3.6) and (3.7) are not 
needed. They are displayed explicitly only for the sake of completeness. 
We now apply this identity to the N terms of (|3.6| ), with 

bf\t)=N{l + a,-'p^,,) (3.8) 



and obtain 

UA{t) 



^ a=l 13=1 1 + a/3 - Jf^a,l3 0/3 " 
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This sum over N terms may be conveniently replaced by a contour integral 
in the complex plane: 

-^w-^/l^(rl^^)^n^ (3.10, 

in which the contour encloses all the a^'s and no other singularity. It is now, 
and only now, possible to let all the a^'s go to zero. We thus obtain a simple 
expression for Uo(t), 

Note that this representation as a contour integral over one single complex 
variable is exact for any finite A^, including = 1. 

3.2 Semi-circle law 

it 

In the large limit, for finite t, the integrand has the limit c'^-") and 
therefore for large N, finite t, Uo{t) approaches 

1 f du it 

By the change of variables 



« = i l-Jl-ll (3.13) 



we have z = ■ Then the integral of (|3.12| ) becomes, after an integration 
by part, 

J Ztxi 



\/-l + -e**^— (3.14) 
V X 27r ^ ' 
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Therefore, we have from ( |3.3| ) 



J Zn 



_1_ / 4- A 
2^V A 

which leads to the expected semi-circle law 

1 



(3.15) 



p(/i) = |/i|pV) = ^\/4-/i2 (3.16) 



3.3 Behavior near the origin 

For A small however, we need to control the large N limit of f/o(t), for 
large t. In this regime the simplified form ( |3.12| ) of ( p.ll| ) is not valid. This 



is why the semi-circle law is modified at the origin. Keeping in mind that we 
shall consider A's of order 1/A^, we write 



P(A^) = iV^ P ^e-^^^Vo(iVV) 

J-oo ZTT 



N 



N dr _ .^^2 1 r du f 1 - 77^ \ 
in which we have defined 

x = NX (3.18) 
In the large A^, small A but finite x, limit, we thus obtain 

N dr r dul .^,^2+ 1 „i 

If we calculate instead 

/(.^) ^ f (3.20) 

we find, after changing u to iu and then r to r — m, 

^ ' J-00 277 J 2TTi ^ ^ 
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Using the standard definitions of Bessel functions, we have 



(-+^) = -^M2x) (3.22) 
^e'-i-H) ^ iJ^(2x) (3.23) 



and thus 



= y^[^'(2a:) + J^(2x)] (3.24) 
Therefore the integration with respect to x'^ is immediate, and we find 

P(A) = ^[Jo'(2iVA) + JK2N\)] (3.25) 

in the cross-over regime in which A^A is finite. The oscillatory behavior we 
expected is described by Bessel functions 

3.4 The edge of the semi-circle 

It is easy to apply this same method for studying the cross-over at the 
other edge of the distribution, in the vicinity of the end point /j, — 2 oi the 
semi-circle. The derivative of the density of state p(/i^) with respect to //^ is 
given by 

r ■. -\ N 

dp{n'^) /■+°° dt f du 



d/i'^ J-oo 27r / 27ri 



-| it 

^ Nu 



e-'*"^ (3.26) 



N{l+u) _ 

Changing t to Nt, and then t to t — iu, and also u to —iu, we obtain the 
factorized expression. 



-iN [ ~ ^^_^)N^-^Nf.^t I rf^.:L^)A^e^^'^'" (3.27) 
J-oo 27r4 -ir J 2nr u ^ ^ ' 



The integration over t is easily expressible with the help of a Laguerre poly- 
nomial: 

n+oo I 



J— oo i — tt 

= -27ri^e-^>''L'j^{Nii^) (3.28) 
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The contour integral over u turns out to be also expressible as a derivative 
of a Laguerre polynomial; we end up with the simple expression for (|3.26| ), 

= -Ne-'^^^L'^iN^j?)]' (3.29) 

Using standard identities for Laguerre polynomials it is easy to verify that 
this leads to 

p{^?) = Ne^'"^'[L^{N^i')L'^^,{Nfi') - L^_^{Nfi')L'^{Nfi')] (3.30) 

The orthogonal polynomial method , which led to the expression (2.8), may 
be cast into this form through Christoffel-Darboux identity. However for our 
purpose, the cross-over distribution near the edge of the semi-circle, it is 
much more convenient to return to the integral In of ( p.28D and to use the 
saddle point method: 

r+oo 

/ dte-^^^fi (3.31) 



oo 



where Sejf is given by 

Seff = -In t + ln(l - it) + i^iH. (3.32) 
The saddle points tc become degenerate at yU = 2, since 



-i ± ./4 - 1 

= ^ . (3.33) 



Then we change variables to 

fx = 2 + A^-"x, 



t = -^ + N~^T (3.34) 



and expand Seff up to r^. This leads to 



7r7 

Soffit) = 2 + - + 2N-^x 
16 

+ i—T^N-^^ + 4:iN-"-^Tx (3.35) 
3 
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We thus find that there is a large A^, finite x limit, provided we fix the two 
unknown parameters a and j3 to 

« = I P=l (3.36) 

We repeat this for the m- integral of ( |3.27| ). We then find that the leading 
terms of ( p.35| ) of order 1, as well as the term 2xN~'^^^, cancel with terms of 
opposite signs in the w-integral. Thus we obtain the following expression for 
the density of state near the critical value = 2, 

= -Nh-|\A,[A^N^^,-2)]\' (3.37) 

where the Airy function Ai{z) is defined as 

AA{3a)~^/^x] = ^ — / cos{at^ + xt)dt. (3.38) 
vr JO 

This Airy function is smoothly decreasing for /i > 2 but it gives oscillations 
for jj, <2. 

3.5 Two-point correlation 

The application of this method to the two-point correlation function is 
straightforward; it leads to a compact and useful integral representation. 
Let us briefiy describe the procedure. Introducing the same source A in the 
probability distribution as in (|3.1|) , we consider the Fourier transform of the 
average two-point correlation: 

f/f (ti,t2) =< lTre^*^^^^lTre^*^^'^ > (3.39) 

from which we shall compute the two-point correlation function p*-^-'(Ai, A2), 
after letting the source A goes to zero. The normalization conditions are 
then 

f/f (tl,t2)=f/f (t2,tl), 

f/f (t,o) = f/J)(t), 

Ull\0) = 1. (3.40) 
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After performing the Itzykson-Zuber integral over the unitary group as in 
(|3.6|), we obtain through the same procedure, 



t/f (ti,t2) = ]^ E / nrfr,|^e-^Sti-(i+-)+^*--+^*-- (3.41) 



ai,a2 i=l 



where we omitted the overall normalization, as we are allowed to do (as 
exoplained earlier). The integration of rj is again done with the use of (|3.42|) , 
in which now 



f,{"i.a2) 



A^(l + a/3 - —{iti6p^ai + it2Sl3,a2))- 



(3.42) 



(2) 

Thus we have the following expression for [/^ (ti, ^2) after restoring the nor- 
malization. 



iV2 

X n 

/3<7 



N N 

E n 

ai,a2=l f3=l 



1 + a/3 



H AT 



a/3 a^ "jY ('^/3ai ^^ai) ~^{^l3a2 '^7«2) 

(a/3 — a^) 



(3.43) 



Keeping track of all the terms in which the Kronecker 6a,p^s do not vanish, 
we obtain 



N 



\\ + aaj(l + a, 



012 , 



N . 



l + a„, -f)(l + a,, 
(a -a _ ih-ii2.) n _ „ - iii)(n -n - ih.) 

\^ai ^a2 jv ' 1 r V ai 7 N /V"'ck2 "7 TV ' 



iV 



iti-it2 \ N 

(c^ai ~ 0,02) 77^(01,02) ^'^"1 ~ a^)(aQ,2 — a^) 



(3.44) 



The Fourier transform of the first term of ( 3.44 ) corresponds to 



1 



dhdt2e-''^^^~'''^'u!>'\h + t2) = -6{Xi - A2)po(Ai) (3.45) 



It could thus be omitted for Ai 7^ A2 but, remarkably enough, the contour 
integral that we shall now consider, will be simpler if we retain this term. 
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Indeed let us consider the integral over two complex variables u and v 

il+u)il + v) 




(l + «-f)(l + i;-f,j 



V — 



N 



(3.46) 



It is straightforward to verify that this expression reduces exactly to (|3.44|) , 
provided we choose the following contours: we first integrate over a contour 
in V which circles around the a^'s and no other singularity. Then taking a 



"02 } 



we integrate over a contour in u which surrounds the 
with «! 7^ these poles yield the sum over ai and 



N ■ 



residue at v 
{N — 1) poles u 

a2 of (|3.44|) . The contour in u has to surround also the pole u 
Remarkably, this last pole reproduces exactly the first term jfU^^ti +^2) of 

We are now again in position to let all the a^'s go to zero: 



1 



tlt2 



dudv 
(27rz)2 

Nv' 



u 



){l+v 



it2 ' 
N ' 



N 



'1- 



Nu 



(3.47) 



The last bracket in ( p. 47] ) is a difference of two terms. Keeping the one in 
this bracket we obtain the disconnected part UQ'\ti)UQ'\t2)) of UQ^\ti,t2). 
The second term gives thus the connected two-point correlation p^^\\i, A2). 

In the large N limit, the connected two-point correlation may be then 
immediately obtained. For finite ti and ^2 we have, in the large limit 

dudv 



UP{h,t2) 

Changing variables to 



N{2m) 



it2 



[u — V] 



_g u(l + ii) v{l-\-v) 



1 



1 



uil + u) 



Z2 



v{l + v) 



and denoting 



u{z) 



(3.48) 



(3.49) 



(3.50) 



14 



we have 



U!:'\h,h) = j^l dzrdz2e-''-'-~''^'Hn[u{z^) - u{z2)]. (3.51) 

Indeed we have used -^-^hi[u — v) = („_^^)2 (^)(^)? ^'^'^ then integrated 
by part over Zi and Z2- We then Fourier transform over ^1,^2^ 

Pf\\iA2) = -]^^35a;1^[«(^i) -"(^2)] (3.52) 

This result could be derived by other methods, and indeed has been obtained 
in a somewhat different form by Ambjorn and Makeenko and others [0 . The 
derivation given here however gives the correlation function even for finite 
A^. One can check easily that the derivation that we have given here through 
Kazakov's representation, could be repeated for calculating the two-point 
correlation function in the unitary ensemble, in which the universality of the 
two- point correlation function has been studied |T2|, |13|, 0, 0, |T^, [T8[. 

It is not difficult to verify that the above expression for the connected 
two-point correlation is a compact version of the one that we would have 
deduced from the orthogonal polynomial method, with the kernel K{Xi, A2) 
in (p.8|) . Indeed, after letting the a-y's go to zero, we determine the connected 
two-point correlation function to be 

, u-iti v-it2 .Ar,l + u.,^,l + v.^ 



1 + u — iti 1 + V — it2 u V 
X — (3.53) 

[U — V — lti){U — V + tt2) 

This expression takes a factorized form if we shift ti to ti — iu and ^2 to 
t2 - iv: 

(Ai,A2) = -I{NXi,NX2)I{NX2,NXi) (3.54) 
in which we have defined 

/(A,.,,) =f^i± f^f fi^fl (3.55, 



27r J 2Txi \t + i 
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Taking the residues at the poles t = —i and f = 0, we find 

N-l 

/(Ai, A2) = e-^^ J2 Ln{\i)Ln{\2) (3.56) 

ra=0 

Thus /(Ai, A2) is simply an integral representation for the kernel K{\i, A2) in 
( |2.8|) . The connected two-point correlation function p^^\n,u) is expressible 
through the Christoffel-Darboux identity, 

- P^ (;,2_^2)2 ^^■^^) 



4 Use of Grassmannian variables 

In the previous section we have used a source representation, which is pow- 
erful and simple. However we would like to investigate the same question of 
the edge behavior near the origin for more complicated ensembles of block 
matrices which arise when the randomness is due to random couplings be- 
tween neighbors on a lattice. Unfortunately we have not found any simple 
extension of Kazakov's method and we have to use Grassmannian variables, 
as often in disordered systems, to solve the problem 0. Before going to a 
lattice of matrices we return to the case that we have solved in the previ- 
ous section, in order to show how to recover the same results from this new 
method. We shall then extend it to lattices in the next section. 
We begin with the identity 

N 

M-f^ = -iN / Y[{ducduldvcdv*) 
1 

X ulu,eM^Y.^KMcdUd + vlM,dVd)] (4.1) 

c,d 

in which the m's are commuting variables and the f 's are Grassmannian. 
Indeed with the normalization 

/ dvdv*vv* = - (4.2) 
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one verifies that 

jdudu*dvdv*exp[—{u*u + v*v)] — 1 (4.3) 
In order to apply this to the matrix 

which is 2N x 2A'", we have to decompose the 2A'"- component vectors u and 
V, into A'"-component vectors; 



(4.5) 



and we obtain 



1 1 i f ^ 

Tr = — / Yl (dacda*.dbcdb*.dacda*.d3cdd*) 

2N z-M 2 7 c c c c yc yci 

[(a* ■ a) + {b* ■ 6)]exp(iiVz[(a* ■ a) + (6* ■ h) + (a* • a) + • 
-z7V[a:C,% + alCl(5a + 6:Ceda, + /3:C^«,]) (4.6) 

The average over the Gaussian distribution for the matrix C 

P{C) = ^exp(-ArTrC^C) (4.7) 

Zj 

is then easily performed, since 

< expiArTr(AC + //C^) >= exp(-A^Tr(AAi)) (4.8) 
This gives for the average resolvent 

G{z) = < -^Tr^— > 

i f ^ 

= -- / l[{dacdaldbcdbldacdaldpcdp*)[{a* ■ a) + {b* ■ b)] 

exp(i7V^[(a* • a) + (6* • b) + {a* ■ a) + (/5* • /5)] - iV[(a* • a) (6* • b) 
-{a* ■ a){p* ■ (3) + {a* ■ a){b* ■ (3) + (/3* • b){a* ■ a)\) (4.9) 
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Note that there is a minus sign in front of the four Fermi interaction due 
to the Grassmannian algebra. This four Fermi term may be replaced by an 
additional integration, since 

TV r 

expAr(a* ■ a){(3* ■ (3) = — d^aexp {-N[(t*(t + a* {a* ■ a) + a{(3* ■ (3)]) 

TT J 

(4.10) 

Substituting this into ( [4.9|) we can now perform the integration over the 
anti-commuting variables, 

N 

/ \{{dacdald(3cd(3l)e-xY>{-N[{a' ■ a) (a* - iz) + {(3* ■ (3){a - iz)] 

•' c=l 

-N[{a* ■ a){b* ■ p) + {p* ■ b){a* ■ a)]) 

= {^rdet(,f-'' i^x.^i 

vr \\o >< a\ a — tz 

= (— )''^[(^* - ^z){a - iz)f-^[{cr* - tz){cr - %z) - (a* ■ a) (6* ■ h)\ 

71 

(4.11) 

We are then led to 

G{z) = - — (— )2^ / d^aT\{daJaldb,dbl)[{a* ■ a) + {b* ■ b)] 
2tt n J ^ 

[{a* - iz){a - iz)f-'[{a* - iz)ia - %z) - [a* ■ a){b* ■ b)] 
exp{iNz[{a* ■ a) + {b* ■ b)] - N[{a* ■ a){b* ■ b) + a* a]) (4.12) 

The integrand is a function of the lengths of the complex vectors a and 6; 
integrating over the angles: 



YlidaJaDfia*, a) = dxx'^-'fix) (4.13) 



we end up with 



[{a* - %z){a - iz)f-%a* - iz){a - iz) - xy\ 
exp(A^[z2:(x + y) - xy - a*a]) (4.14) 
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This representation of the average resolvent in terms of an integral over four 
variables is exact for any N. 

Let us note that, if we had calculated, instead of the average resolvent, 
the average value of 1, we would have obtained by the same method, the 
identity: 

' dPa dx dy{xy)^-\a* -iz){a -iz)]^-^ 



TT [N\f J Jo Jo 
[{a* — iz){a — iz) — xy\exjp{N[iz{x + y) — xy — o"*o"]) = 1 (4.15) 

This identity should manifestly hold, but it is not quite trivial to check it; 
since it provides an interesting verification of the consistency of the formation, 
let us note that one can derive it by replacing, in the l.h.s. of ( [4. 151) , the 

bracket [(cr* — iz){a — iz) — xy] by 

1 92 d"^ iz d d d d (A^a\ 

followed by integrations by part and a few lengthy manipulations. This 
identity will be useful in a few moments. 

The expression (|4.14| ) of G{z) is well suited to study the large A^-limit. 
The integrand involves a factor exp(— A^S"), with 

S{a, a* , X, y) = —iz{x + y) + xy + a* a — \n.[xy{a* — iz){a — iz)] (4.17) 

and the large N limit is therefore governed by the saddle-point at which 5* is 
stationary. The equations for the saddle-point lead, away from the vicinity 
of 2; = 0, to the equations 

xc = yc = (yc = < = \[^z + v/4^] (4.18) 

The sign of the square root is chosen so that the imaginary part of G{z) 
above the cut is negative; note that in the saddle-point method a and a* 
become independent complex variables. In order to obtain G{z) with the 
proper normalization we should include the Gaussian fluctuations around 
this saddle-point. It presents no difficulty; however we can bypass the whole 
calculation, if we note that the integrand in ([4.14|) differs from that of ( [4.15| ) 



by a factor —i{x + y)/2; if we had calculated ( |4.15| ) by the same saddle-point 
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technique we would have obtained one, a very cumbersome method to get 
one indeed. But this immediately tells us that, for N large. 



G(z) = -'-(xc + Vc) = - iV4^] (4.19) 
from which one recovers the semi-circle law 

p(A) = --ImG(A + iO) = — \/4^^ (4.20) 

TT 2tT 

However this solution is not valid near the vicinity of the origin. Indeed 
the saddle-point equations gives z{xc — Vc) = z{af. — a*) = 0; another way of 
realizing that there is a problem near 2; = 0, is to calculate the determinant of 
the Gaussian fluctuations near the saddle-point, which vanishes when = 
1, i.e. for z — ±2 or z — 0. Near the end-points of the semi-circle the 
phenomenon is well-known and leads to an Airy type cross-over function in 
a region of size A^~^/^ near the edge, which we have discussed in the previous 
section. Near the origin, we need a separate analysis. We will focus on a 
range of size 1/iV, in which the variable 

C^Nz (4.21) 

is finite. We write a — u + iv, a* — u — iv and translate the real u- contour 
by iz; this gives 



2iAr-i 



GiC) ^ -^ j^ J dudv dxj^ dy{xyr-\x-ry)[u'-rv'\ 

[u^ -\-v^ — xy]e'xj^{—N[xy v^])exp(iC(x -|- y))e'x^{—2iuC, + ^C^) 

(4.22) 

Going into radial variables for u and and changing y to p and q with 
xy — p, X + y — 2y/pq (Jacobian J — {(f — 1)~^^^, domain of integration 
p > 0, g > 1), we obtain easily, if we drop terms of order 1/N, 

exp(— iV[p + r — 2 — ln(pr)])exp(2i^(7i/p)exp[— 2iCA/rcos6'] (4.23) 
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The remaining integrals over p and r are governed, in the large N limit, by 
their saddle-points at p = r = 1. However the presence of the odd factor 
(r — p) in the integrand implies to expand beyond the leading non- Gaussian 
order. Changing variable q to cosh0, we integrate and 6, 

POO 

c/0cosh0e2<v^'=°^^^<^ = i^i(-2<VP) 

= -|[Ji(2iV/^v^)+ziVi(2iV/i^)] (4.24) 

in which N ^ is the real part for i.e. the eigenvalue of the matrix multiplied 
by N- 

I ^exp(-2iA^/iv^cos^)de = 2^Jq{2N ^ly/^) (4.25) 
Jo 

we have for the imaginary part of G{Q 




exp(-A^[p + r - 2 - \n{j)r)])Ji{2N ^i^)Jq{2N ^i^) 

(4.26) 

Expanding the Bessel functions Ji{2N Hy/p) and Jo{2Nfiy/r), with r = 1 + r' 
and p = 1 + p' up to linear order, we get 

Ji{2N^i^) = Ji(2iV/i) +p'[N^iJo{2Nfi) - ^Ji(2iV/i)] (4.27) 

Jo(2A^^v^) = Jo(2Ar/i) - N^r'Ji{2N^) (4.28) 
Then, by the Gaussian integration over p'^ and r'^, we have 

P(/^) = ^[Jl{2N^^) + Jl{2N^,)] (4.29) 

Since Jq{x) ~ 1 and Ji(x) ~ x/2 near x = 0, the density of state p(/i) is 
proportional to /i for small /i. The density of states that we have found 
agrees with (|3.25|) . If we average over these oscillations with the appropriate 
width, we find < p{^) >= I/tt, the value that it takes at the origin in the 
saddle-point method for the large limit ([4.20|) . 
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5 Ring of matrices 



We now extend the previous Grassmannian method to the case of a lattice 
of matrices. We consider a ring of L points, with L even, in which the 
neighboring sites are coupled by complex N x N matrices. The previous 
section corresponded to L = 2, a lattice of two points, coupled by the matrix 
C for one orientation of the link, and for the opposite orientation. The 
simplest extension, L = 4 consists of a AN x AN matrix given by 



M 



( ° 


c\ 










c\ 




















c\ 

/ 



(5.1) 



where Cj is a x complex matrix. This matrix M represents a random 
hopping between the nearest neighbour sites of a lattice of four points on 
a ring. The disorder is off-diagonal, and the hopping terms are N x N 
random complex matrices. This matrix M has also pairs of opposite real 
eigenvalues, since it anticommutes with the "75" matrix. However this L = A 



case (5J.) happens to be reducible to the previous case through the orthogonal 
transformation which exchanges indices two and three; namely if we change 
M to P-^MP with P(l, 1) = P(2, 3) = P(3, 2) = P(4, 4) = 1, and all other 
P{i, j)'s equal zero, one sees easily that the problem is mapped into the L = 2 
case with replaced by 2A^. This is not true however for larger rings. 

For L even the matrix M has again pairs of opposite real eigenvalues since 
it anti- commutes with the "75" matrix made of L block matrices consisting 
successively of the unit matrix and of minus the unit matrix. One can thus 
consider again the problem of the behaviour of the density of eigenvalues in 
the scaling range near the origin. We shall prove that the previous result 
still holds up to scale factors. 

The formulation of the previous section for 2A^ x 2A^ matrices, may be 
easily extended to the LN x LN matrices, corresponding to a one dimensional 
lattice. The 2A^ component vectors u and v in ( [4.5|) are now LN- component 
vectors; they are conveniently decomposed into iV-component vectors, ai,ai 
{i = 1. ■ ■ ■ , L) where the a's are Grassmannian variables: 



u 



(5.2) 
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Since the matrices Ci{i = 1, ■ ■ 
the sites i and i + 1 , we have 



L) are associated with the hopping between 



1 



-Tr- 



1 



LN z-M 



^ i=i 



li ■ ■ ■ OtL 



[{a\ ■ ai) + ■ ■ ■ + (a^ ■ aL)]exp(iA^z[(a* ■ ai) + ■ ■ ■ - 
+ {al ■ ai) H h (a2 ■ (^l)] - iN[a\^{Cl)cda2d + 



1% ■ OlJ 

• + al^{CL)cdOiid]) 
(5.3) 

The average over the Gaussian distribution for the matrices Ci, ■ ■ ■ , C^; 



1 



P(C) = ^exp(-iVTr[CjCi + ■ ■ ■ + ClC^]) (5.4) 

leads to the average Green function G{z). As in ( ^.9| ), we replace the four 
Fermi term by an integration over additional variables cr^'s, and, as in ( d.lOj) , 
we integrate over the anti-commuting variables: 

T = y daida\ ■ ■ ■ (iaLda2exp(— A^[(q;^ ■ ai)(a^ ai — iz) + ■ ■ ■ 
+ {al ■ aL)(a2 + <tl-i - iz)] 

-N[{al ■ ai){al ■ ^2) + ■ ■ ■ + (a^ ■ aL){al ■ ai) + C.C.]) 



det 



/ Si 

02 >< o-il 



X 



where 



\ \aL >< ai| 



\ai >< 02! 












IttL >< CtL-ll 



\ai >< 


St, 



TT 



cilI \ 

/ 

(5.5) 



5i 
^2 



5/ 



(5.6) 



In the calculation of the determinant, the vector appears only through its 



squared norm denoted as Xi 



a,- 



Then the Green function G{z) becomes 



L TT 



Y[ d^adaida* (^ Xi)T[a, 



X) 



i=l 
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exp{iNz ^ — A^[^ cr*(Tj + X1X2 + ■ ■ ■ xlXi]) (5.7) 



where the factor {N/tc)^ is due to the introduction of the o"j's (^4.1CI|) . Chang- 



ing variables from to (|ajp = Xi), as was done in we have 

(xi ■ ■ ■ xl)'^~^{xi H h XL)T{a, x)exp[iNz{xi H h x^) 

-A^^a*(Ti-A^(xiX2 + --- + XLXi)] (5.8) 

In the large N limit, repeating again the argument of the previous section, 
we have a saddle-point 



/ \ r *\ r \ iz + ,_ 

(^i)c = = (ai), = (5.9) 



and it leads to a semi-circle law for the density of state (note that ( |5.9| ) differs 
from (|4.18| ) since we have used in this section a slightly different normalization 
of the probability distribution). 

As before the determinant T{ac,Xc) vanishes at 2; = . We then have to 
expand the variables of integration around these saddle points. We shift the 
complex variables a* — iz/2 to a* and cij — iz/2 to (Xj. This gives 

Si---Sl = (a* + az.) ■ ■ • {al + a^^i) (5.10) 

It is convenient to replace the variables Xj's, which are positive by definition, 
to the variables (Ai, A2, ts, t^,---, ti) defined by (i = 1, ■ ■ ■ , L), 

Ai = Xi+XsH hXL_i 

A2 = X2+Xi^ Vxl (5-11) 



and 



Xi = Ai(l -t-i tL^i) 

X2 = A2(l - t4 ^l) 

X2n-1 = Alt2n-1 (n ^ 1) 

X2n = A2t2n (^ 7^ 1) (5.12) 
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The Jacobian J for this transformation is J = (AiA2)^'^^ ^. The Green 
function is then written as 

G{z) ^ -- j n^^V, j dXid\2dh---dte{XiX2)^-' 

[AiA2]^[(l - ts tL-i)(l -U tL)t3U ■ ■■tLf-'iXi + A2)T 

exp{iNz{Xi + A2) - Ar(AiA2)/i(t)) 
^ zN ^ 

exp(-7V ^ a*a, - i— ^{a* + a,)) (5.13) 
i=i ^ i=i 

in which 

h{t) = htL + t2ti + + ■ ■ ■ + tLtL-1 (5.14) 

^1 = 1 — — t5 — • • • — ^L-i 

^2 = 1 - ^4 - ^6 tL (5.15) 

We return now to the variables p and q as before, 

A1A2 = p 

A1 + A2 - 2v^g, (5.16) 



with Jacobian J = l/\/g^ — 1. In the large N limit, the saddle point for the 
tiS is 

2 

^2n-l — 'l^2n = (5-17) 

and h{tc) = Since T{ac,Xc) is vanishing, we have to expand ai and 

p, g, ti around their values at the saddle-point. First, we integrate over q 
with the change of variable q — cosh0. Then we get, 

^^^^ ^ ~L~(my~ J ^ J ^yn^^^ 

[(1-^3-^5 tL-l){l -h-te tL)t3t4h ■ --tLr-^T 

i^i(-2zCVP) exp[-Nphit) - iV^a*a, - ^CE(^i + ^«)] (5-18) 

The saddle-point is now tc — Pc — ^ and h{tc) — j. The determinant T 
vanishes at this saddle-point. With the polar coordinate representation for 

CTi = ^ie'^' (5.19) 
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the saddle-point appears now at 

{02n-l)c = -{02n)c = ^1 (5.20) 

We denote the deviations from this saddle-point as 9'^, 

02n~l = 01 + ^2n-l 

= + 9'^^ (5.21) 

with 9[ = 0. The determinant T is independent of 9i, and thus the integral 
over 6*1 which appears in ( ^.18| ) is 

A = J d9ie-"^^ti v^.cos{e,-{-iye[) ^5 22) 

Expanding up to second order in the ^^'s , after integration over ^1, we get 

A=.27rJo(CEv^)+0(^^f) (5.23) 

1=1 



It is only the first term of ( p.23| ) which matters, since T is vanishing for the 
saddle-point, and the fluctuations over the ^^'s are of higher order. Therefore 
the imaginary part of G{z) is given by the product of Jq and Ji as before. 

ImG(.) = — (^,). il^f j "gdnWde', j dpp- 



i=l j=2 



L 



J n dU{ht2 ■ ■ ■ tLf-^6{ti + ta + ■ ■ ■ + ti-i) 

i=l 

6{t2 + h + --- + tL)T 

exp(-iVp/i(t) -iVX:r,)Ji(2Cv^)Jo(CEv^) (5.24) 



i=l i=l 



The factor (2) comes from Jacobian of dr. 



We still have to integrate over p, ti,9[,ri. However, the integration around 
the saddle-point for ti and 9' may be avoided, although they could also be 
obtained after quite lengthy calculations. Indeed one may notice that if we 
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had calculated, instead of the average resolvent, the average value of 1, we 
would have had an integral similar to (4.15). If we attempted to compute 
one through a large A'"-limit analysis, we would obtain again a product of 
two Bessel functions Jo{NnJ2 \ff\) JiC^^l^y/p), in which Nfi is the real part 
of (. The coefficient of this term should thus vanish. Therefore if we replace 
p and Ti by their saddle-point values, then the contributions of the fluctua- 
tions of 9[ and tj, which also have a factor Jq[LN(i/ \/2)Ji{LNfi/ \/2), should 
cancel each other. The only difference between the calculation of the average 
resolvent and of the identity, is a relative factor ^/p. This contribution is also 
exactly cancelled by those coming from the expansion of the Bessel function 

Ji(2Ar;,^) ^ Ji(^^/.) +p'^^Jo(^) -p'_J,(^) (5.25) 

It is thus possible to replace ti and Qi by their values at the saddle-point. 
The Bessel function Jo(C 5] y^) expanded as 

J„(iV;,i:v^J)-^o(^)-^i;.:i.(^) (5.26) 

We are thus left with two terms, which are proportional to \Jq{LN 
and \J \{IjN \x I . The coefficients of these therms are evaluated by inte- 
grating over p' and r' . The determinant T factorizes as 

z. 

It is simply the product of the eigenvalues for a periodic chain with nearest 
neighbor interactions. Setting p— ^ and expanding this product up to 
order p', we find 

V(;^)2^ n (sin^(|^) - y^o.^(^^k)) 
--P'(;^)2^ nsin^(^A;) 

= -pll?-^ (5.28) 
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where we have used the identity: 



n ^H-k) = — T (5.29) 

k=i "' 



The o integral, which appears in this calculation, may be done exactly. After 
integration over the crj's, we have 

L 

/tvl = y n ^'^^[(^1* + ^^^)(^2 + fxi) ■ • • {ol + (Ti-i)]^ exp(-iV^ 



TT 



^C(iV,L)(Ar!)2 (5.30) 



The quantity C {N, 2k) is expressible as an integral over angles; for instance 
in the k = 3 case, it reads 

C(iV, 6) = / / de,de2[l + e'''f[l + e-^'^-^c^ll + e'^'f (5.31) 

(/vrj^ JO Jo 

In the large N limit, we exponentiate the integrand of ( p.31| ) and expand the 
^'s around 9 = 0. Then we have 2(8^)/{tt\/3N). For general k, we get 

C(N,2k) - , , (t7)^ , — (5.32) 



From the identity. 



ncos(-) = ^ (5.33) 



and setting = ^, we have 



2 • 



n(i+cos(^))=2#-n^o«'(^) 

n=l ^ n=l ^ 

L 



Thus we get 



2-2L (5.34) 



2^ 2T 1 

C(iV,L)^^^^-^ (5.35) 
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and Inl in ( |5.30| ) reads, for the case of an arbitrary L , 



i-NL 



(5.36) 



The remaining integration over p may again be done by the saddle point 
method; setting p = 4- + p', we have 



Ip = J dpp 2 2p'^e ~ (— ) 2 2e 2 y dp'p'^e 

VZTT .Zv^, NL 7 LN , 

Thus the piece of the imaginary part of G, which is proportional to Joi^^C)^ 
is obtained by mutiplying together the coefficients of ( ^■25|) , (|5^28|) , (|5.36|) ,( ^.37| ) 
and the contribution coming from the integrals over the tj's , which is given 
in an appendix: 

(ImG')^. ^ --{^)DnNf,j',{^Nf,) 

= -^nN^^J',i^N^^) (5.38) 



where Nfi is a real part of (; the calculation of D is given in Appendix A, in 
which it is shown that D = /L^. We have also the same coefficient for 
the other piece of the imaginary part which is proportional to Ji{LNn\/2/2). 
Thus we obtain for this one-dimensional ring, 

= -/^[Joi^^t^) + JK^^f^)] (5.39) 

If we replace LN ^/ \/2hy 2x, then this equation becomes p(/i) = -^^[Jq{2x) + 

Jf{2x)]; it is identical to ( ^.25| ) (except for a trivial factor 1/a/2 which comes 
from a different normalization in the probability distribution in this section; 
for that reason the edge of the density of state is now at Zc = 2^2 instead of 
2). Thus we have obtained a behavior at the origin for a chain of L matrices 
which, up to a normalization, is identical to the previous simple case. 
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6 Lattice of matrices 



If we now consider a higher dimensional lattice, with a total number of lattice 
points equal to L, and periodic boundary conditions, one may ask again the 
same question. As before, we take an LN x LN random matrix , with N x N 
block elements, corresponding to the hopping between nearest neighbours on 
a lattice. Non-neighbouring sites are not coupled and are represented in the 
total matrix by block matrices of zeroes. Generalizing the expression of (|5.3|) , 
we have 



LN 



J- If T-r 

-Tr— ^ = -- y J] dau ■ ■ ■ dal^dau ■ ■ ■ da^^ 

[{a[ ■ ai) H h (a^ ■ aL)]exp(iiVz[(a^ ■ ai) H h (a^ ' 

+ {al .ai) + --- + {al- az.)] - iN[al{Cya2 + •■■]) (6.1) 



The last term of (6J) reproduces the connectivity of the lattice. 

Integrating over the random complex matrix Cjj, and using the vari- 
ables of ( [4.10| ), we have, as in ( |5.8| ), 



<^,3> 



(xi ■ • ■ xl) (a;i H h XL)T{a, x)exp[iNz{xi H h xl) 

E - N{xMjXj)] (6.2) 

<i,j> 

where we used the notation |ajp = Xj and < i,j > are a pair of nearest 
neighbours; Mjj is the adjacency matrix of the lattice. One can use a method 
similar to that of the one dimensional chain; we divide the lattice points into 
two groups, with odd and even indices. This change of variable leads to an 
expression similar to ( p.4| ) 

i A^2L^L(Af-l) r r l_ 

G{z) = — J d^ai J d\idX2dti---dtL{\i\2)^ ^ 

i\l\2)^[tlhts ■ ■ ■ tLf-\\l + \2)T 

S{h + ts + ■ ■ ■ + tL-l - l)S{t2 + U + ---+tL-l) 

exp{tNz{X, + X2) - NX,X2h{t) - N J2 (6-3) 



'■,3 

<i,j> 
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where h(t) = UMijtj. By the change of variables, A1A2 = p and Ai + A2 = 
2y/pg, we integrate over q, and obtain immediately Ki{—2iC,^) as (2]2^). 
We replace Nz by C- The saddle-point is (ti)c = |;, Ac = f, h{tc) = f and 
Pc = Y^- (Note a factor of two in the normalization compared to the one- 
dimensional problem). We shift the complex variable (J*j — iz/4 to a*j and 
(Tij — iz/4: to CTjj. This shift gives an extra term to ( |6.3| ) of exp{—izNJ2{<^lj + 
(Tjj)/4). We write cxjj as 

^ij = V^e"'"'' (6-4) 
We expand the variable 9 around the saddle point 9ij, 

02i-l,j = 6*1,2 + ^2i-l,j 

02i, = -01,2 + 9',,^^ (6.5) 

As in the previous cases, we may simplify 

exp(--C5]Ki + -exp(--C5Iv^cos0i,2) (6.6) 

The deviations from the saddle-point ±^12 may be neglected since T is van- 
ishing at the saddle-point. Therefore, by integration over ^12, we obtain 
again Jo{^J2^/f\j), where Nfi is the real part of (. We note the values 

at the saddle-point are {rij)c = \ and pc = Thus we have the following 
expression: 

27r2 N'^L L{N^i) I . r NL 1 f 

ImG = q)'J Udr.,UdO:J dpp^ J YldU 

[ti ■ ■ ■ tLf-'T5iti + ts + • ■ ■ + tL-i - l)5(t2 + t4 + • ■ ■ + - 1) 

,iV/i 



M2Nf,^)JoC-^ E v^)e-^^'^(*)-^S-'^ (6.7) 

<i,j> 

The determinant T is expanded around the saddle point. If, for defi- 
niteness we specialize to a two-dimensional square lattice, T contains the 
factor: 

n [1 ~ yVPCos^ki - yv^cos-^/ca] 
fci,fc2=i ^ y/L L ^ L 

4 -n- . 2 ^ 27r , 2 ^ 27r , , , , 

= (1 - y v^) 11 [1 - y v^cos^/ci - -^cos^fcaj (6.8) 
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Expanding p around the saddle point, p = ^ + p', we find that the first 
factor of ( |6.8|) becomes 

1-^Vp=-^p' + 0{p") (6.9) 
We also have to expand the Bessel functions and we obtain 

M2Nfi^) = + -j^p M^) - J2P M^) (6.10) 

The last term of order p' cancels as before with the term of l/i/p in (|6.7|) , 
as seen from the expansion 

VP^^i^ + ^P') (6-11) 

Repeating the procedure used for the one- dimensional chain, we obtain an 
identical form for the density of state near the origin, 

p(/.) = C/i(Jo^(^^/.) + Ji^(^/i)) (6.12) 

where the coefficient C is For this two-dimensional lattice, the end point 
of the semi-circle density of state is, with our normalizations, /i = 4; this is 
why we have an extra factor of one half in our scaling variable compared to 

7 Non-Gaussian probability distribution 



Up to now, we have considered a Gaussian distribution for the random matrix 
C ( ^.21 ). If one modifies this distribution, for instance by adding quartic 
terms in the exponential, the average density of eigenvalues will no longer 
obey a semi-circle law. However earlier studies on the usual unitary ensemble, 
revealed that the behavior near the edge of the semi-circle was not affected 
by non-Gaussian terms. The cross-over there is always given in terms of Airy 
functions ( p.37| ). We are thus led to investigate whether this universality is 
also valid for the oscillation of the density of state near the origin that we 
have found for the block-matrices that we are considering in this article. 
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Let us consider a non-Gaussian distribution, for the simple one-matrix 
(i.e. L = 2-model) of section 3, but it will be clear that the proof of univer- 
sality will apply to any lattice of matrices. Consider for instance the following 
P{C): 

P{C) = \ exp(-ArTrC^C - N gTY{C^ CC^ C)). (7.1) 

Zj 

Note that the factors of in the exponential are such that the average 
correlation functions of the eigenvalues, and in particular the density, have 
a finite limit when goes to infinity. Had we put a higher power of in 
front of the quartic term, then it would not be the case; with a lower power 
of A^, it would not contribute at all. Then in the cross-over region of size 
near the origin, one could argue at first sight that the quartic terms modify 
simply the overall scale, but not the oscillatory behavior. Indeed let us write 
the average density 

Po(r) =< ^Tr5(r - C^C) > (7.2) 

We then integrate over the unitary group and scale r, as well as the eigen- 
values of C^C by . Then one sees immediately that the quartic terms of 
( |7. 1|) give correction of order . However the overall normahsation remains 
^f-dependent and therefore it modifies the scale of the cross-over function by 
a factor which is the ratio of the partition functions for the g ^ problem , 
and the Gaussian one. 

But this simple-minded analysis is based on letting A^ go to infinity first, 
and in fact it is slightly misleading. As will be seen now, the question is 
more subtle, and the non-Gaussian terms affect more than simply the overall 
normalization of the density, as we pretended in the previous argument. 
It has already been found in similar problems [0, that by letting A^ go 



to infinity first, one computes only a smoothed average of the correlation 
function and this is not what we are considering. Indeed if we smoothed out 
the oscillatory part of the density near the origin, the simple universality 
claimed above would be true : the non-Gaussian part would change only 
the normalization. However, since we are interested in those oscillations, the 
previous argument is not sufficient, and we shall argue now that the non- 
Gaussian terms do modify the period of these oscillations. This change of 
the approximate period of oscillations, is in fact expected; indeed there are 
A^ eigenvalues distributed between zero and the endpoint. There are thus 
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oscillations in the density. If the normalization is changed, and say the 
value of p(r) at r = multiplied by a factor c, then the approximate period 
of oscillations has to be divided by c. 

This may be checked explicitly by returning again to the formulation in 
terms of contour integrals, which we have developed in the section 3; we shall 
apply it now perturbativcly for this non-Gaussian distribution. 

For the non-Gaussian distribution, UA{t) in (3.2) is given by 



UA{t,g) 



/ e-^E('^^+'^««^+s'-|) n,<,-(n - Tj) n dr, 



This is expressible as 



where 



n9 a„ 



7f 



(7.3) 



(7.4) 



N 



a=l 



D^F{t^ 0) 
Expanding this UA{t, g) up to order we get 



(7.5) 



9 



N 



UA{t,g)c^UA{t,g^O)-j-Y: 

-'^ i=l L 



,d^F^ 1 ,d^D^ F 



Noting that 

OUa 
dai 



+ 0(/) (7.6) 



^2 



dai 



;UA{t,0) 



dF 1 _ F dD 

(^)- - -(^) + 2^(^f - 2(—)(—)- 



dF^,dD^ 1 

(7.7) 



the term of order g, denoted by 5Ua-i becomes 



N 



d 



2g ^.dUA.,d\YiD 



1=1 



da. 



N dai 



)(- 



dai 



-) 



(7.8) 
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The expression of D is given in (3.7) and we have 

dl^^y^ ^ (7.9) 

dtti f^. {ai -Qk) 1 + ai 

Using the contour representation in (3.10), we get 
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& 



Noting that 



UAit,0) =-|/^(TS%)'^(^n,^.(^) (7.10) 



1 it , 1 it 



{u - aj)2(ai - ttk) N{u - ak) {u - akf{ak - a^) N{u - a,) 

it 

^ (7.11) 



2u — tti — ttk 

\2(,, _ n, \2 



we are able to express the contribution 6Ua in the contour integral Thus 
we get the contour representation in the first order of g; by letting goes to 
zero, 



2g f du ( 1 + u \ \ I u 



,'^fdfu_( 1 + u \^ iV-1 7,-1 _^ 

it 



U 
N 



u 



Immediate checks for A^=l, 2 and 3, and also the large N hmit, are easy and 
they agree with the result known by other methods. Thus this expression is 
exact to this order in g. 

We now consider the cross-over region near the origin; replacing t by N'^t, 
and u by Nu, we have 

N-1 1 2u- it 1 1 

^ 2m ^ {u - it)^ ~ N u{u - ity ^ ■ 
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Only the last term in the bracket contributes in the large N limit. It gives, 
by the Fourier transform in (3.3), 



5p(/i) = 2gNJl(2x) (7.14) 

where x = Nfi. The overall factor due to the change of normalization is 
1 + 2g, which agrees with the expression for the density of state which was 
found in [2]. Then, up to this order, we may interprete this result as reading 

Po(/i) = Cig)^[.P,{2Nf,^)) + Jl{2N^,^))] (7.15) 

where C{g) = 1 + 2g + Olg"^). This expression makes it clear that the 
integrated density of state remains properly normalized to one. Indeed, ex- 
panding this expression for g small, we find the term gJl{2x) cancel, and 



obtain ( |7. 14| ) . This result is expected; it is exactly the universality which 
was claimed earlier up to this order in g. 

We are thus tempted to conjecture that these oscillations are indeed uni- 
versal, namely that for an arbitrary probability distribution of the form 

Zj 

with C defined over a lattice the density of states near the origin is given by 
p{^^) = —^[Jl{2NLF^,) + Jl{2NLF^i)] (7.17) 
with F{y) some functional of V . 



8 Oscillation and cross-over behavior for real 
matrices 



It has long been known that for real symmetric matrices the relevant Jacobian 
involves the absolute value of the van der Monde determinant and thus the 
corresponding orthogonal polynomial analysis becomes quite complicated. 
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However, by using some remarkable identities, we can actually treat the cross- 
over behavior near the origin of the density of state of a {2N + 1) x {2N + 1) 
real matrix M, which is made of a rectangular {N + 1) x N matrix C: 

M=(° (8.1) 

where is the transpose of the matrix C. As we will see, the cross-over near 
the origin shows a different universal behavior from what we have studied 
earlier. 

To understand why we treat the case of C being {N + 1) x N and to get 
oriented, let us do a simple exercise in power counting. Consider an M x 
real matrix C with its MN real variables. Denote the eigenvalues of the 
N X N real symmetric matrix C^C by r^. We would like to have 

AT 

~ Jl In - rj\dri ■■■drN (8.2) 

To see if this is possible, let us do dimensional analysis and count powers. The 
left hand side has the dimension of r 2 ~ qN'^+n.^ ^^^q other hand, 
the right hand side has the dimension of . Equating MN = N"^ + N 
we find M — N + 1. This explains why we chose C to be (iV -|- 1) x A". 
We can of course also check this Jacobian by direct computation, using the 
Fadeev-Popov method for example. 

Note that if we had chosen C to be A" x A", we would have 

N N 

dC ~ Yi, ~~i=^ n l'^* ~ rj\dr\ ■ ■ ■ drN (8.3) 

i V^i i<j 

The presence of the square root factors can be deduced by dimensional argu- 
ments or determined by a direct computation of the Jacobian. These square 
roots make the calculation much more complicated. See below. 
We have as usual the density of state 

and from the block structure of M, 

pi^i) =< ^Tr5i^i' - C^C) > (8.5) 
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(Note that the {2N + 1) x (2iV + 1) real matrix M we started out with 
has an eigenvalue at zero; obviously, the corresponding eigenvector is the 
vector orthogonal to the columns in the x + 1 matrix C'^. This 
eigenvalue leads to an additonal delta function at the center of the spectrum 
in the density of states.) As in the complex matrix case, we define p(A) =< 
;^Tr5(A — C'^C) >, with p(/i) = yup(/i^). With the probability distribution 



P[C) = -expi-NTiC^C) 

Zj 



(8.6) 



this model is known as the orthogonal Laguerre ensemble pO . 

In an obvious extension of Kazakov's method we introduce an external 
source matrix A, which we take to be an A^x A^ Hermitian matrix diagonalized 
by the unitary matrix U . The probability distribution is 



Pa{C) 



Za 



exp(-A^Tr(:7^C - NHxAC^C) 



(8.7) 



Since C^C is a real symmetric matrix, it is diagonalized by an orthogonal 
matrix O, and then 



TrAC^C = Trf/" 



/ai 



V 



/ri 



V 



O (8.8) 



Since there is no known analog of the Itzykson-Zuber identity for inte- 
grating over orthogonal matrices, we would have been stuck at this point. 
The crucial observation is to notice that we could integrate over the unitary 
matrix V = UO^. Strictly speaking, A is not an external source since we 
integrate over all matrices unitarily equivalent to A. However, since we set 
the Oj's to zero at the end, this procedure gives us the correct value of Uo{t). 
Thus, integrating 



„ „ AT 1 TV 

UAit) ^ dV dr^Yl |ri - rj\{- J2 e''-")e-^S^'-^T^^"'^'^-«^^^^*'" 

i<j a=l 



(8.9) 
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using the Itzykson-Zuber identity, we obtain 

TV 



N TV 

X exp(-iV^ri -iV^r^aj). (8.10) 

i=l 3=1 



where UAit) is normahzed as Ua^^) = 1- 
A remarkable identity 



TV TV 



/ • • • / n ^^^i n signvi - vje 



- J2 ^i^i 



^ 'b,-b 



allows us to integrate over the n's. Let us sketch a proof of this identity. 
The integration region in (8.11) can be divided into N\ regions inside each 
of which the r's are ordered. Thus, the integral above is equal to 

r+oo r+oo „ 

/ • • • / n dri9{n <r2<rs< < rN)e-^'''^ (8.12) 

Jo Jo 

plus A^! — 1 similar integrals with the rj's permuted and with a suitable 
sign given by Ylii^j signvi — rj. (We will not be concerned with the overall 
multiplicative factor in what follows since it is irrelevant to our calculation.) 
To get oriented, let us do the N — 3 case. Change variables by ri — Xi,r2 — 

ri+X2, rs = r2+X3, to take care of the ordering. We can then immediately 

do the integral above to obtain 

^ (8.13) 



b3{b3 + b2){b3 + b2 + bi) 



We now add the 5 other terms and collect denominators; the sum evidently 
has the form 

(8.14) 



&1&2&3(&1 + &2)(&2 + &3)(&3 + bi){bi +b2 + 63) 
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with some numerator f{b) which we can now determine by general arguments. 
Since our integral vanishes whenever any two of the 6j's are equal, we must 
have f{b) = A{b)P{b) where A{b) = Uf<jibi - bj) is the usual van der 
Monde determinant and where the polynomial P{b) must be symmetric and 
by dimension coutning must be of degree . Thus, P{b) is uniquely determined 
to be 61 + 62 + ^3- We have proved the identity for N = 3. We then proceed 
by induction. Assume that the identity has been proved for some A^. Doing 
the integral for A^ + 1 along the line described above we encounter after the 
first step 

1 

bi{bi + b2){bi + 62 + b3)....{bi + 62 + ••• + &7v)(&i + &2 + ••• + + &jv+i) 

(8.15) 

plus permutations. We now first add the terms obtained by permuting 
61, 62, ,bN, holding b^+i fixed. Then by the inductive process we obtain 

1 ^Njb) 1 

n^i bi ni<j<Nibi + bj) {bi + b2 + ... + bN + bN+i) 

plus terms obtained by interchanging &Ar+i with one of the fej's. The subscript 
on Atv indicates that it is the van der Monde determinant for the first A^ 
6j's. Collecting common denominators and reasoning as before, we find that 
the sum is equal to 

UtVbiUi<j<N+i{h + bj){bi + b2 + ... + bN + bN+i) 

By symmetry and by dimensional analysis the symmetric polynomial P must 
be of degree 1 and hence equal to (61 + 62 + ••• + &Af + ^Af+i)- We have thus 
proved the identity. Note that if a square root factor were present, as in ( |8.3|) , 
we would not have this identity and thus would have difficulty proceeding 
farther. 

Since bi is given by 

k = N{l + a,-^6^,,) (8.18) 

we have 

uaH) = ^ e( . ) n ( , "r~l\~K )i ^ r^r^ ^ ^^-''^ 

iV ^—i -L ~r Uq. -yj^a 7 A'' 7 
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This is expressible in a contour representation: 



UA{t) = --7f TT-iTTT. it) IKt^-— — W 



it J 27ri 1 + u- § 2 + U + ^ 



X 



,2 + 2m- f JI u-ay - f 



n i \ \ n (8.20) 



2 + 2u y u — ay 



Letting all a^'s go to zero, we have 

2 -I- 2?/ — — II — — 

x CXto (8.21) 
2 + 2u u 

where the contour is chosen around u = 0. We are able to check this formula 
by calculating the Fourier transform, 

_|^e-^*^t/o(t). (8.22) 

and verifying that we obtain the same result for = 2 and 3 as we would 
have obtained by directly integrating 

p(ri) = j W_\ri- rj\eyi.-^{-N^ri)dr2- ■ -dr^ (8.23) 

We obtain the semi-circle law just as in the complex matrix case. 

Next, we study the cross-over behavior near the center of the spectrum. 
Changing variables u — >■ Nu and t — > N'^t we obtain in the large N limit, 

/ N N f du 2 2 ^ u . , 

Uo(t)- (p e-"— 1 + 8.24 

2itJ 2-ni ^ u - it' ^ ' 

By following the similar procedure as ( p.20|) and ( p.21| ). We find 

-^p(x) = -\ji{xf + -J,{x)Ux) (8.25) 



where x = \/2NX. Thus we obtain 

p(A) = Jo{V2N\Y + Ji{V2NXY - -yl — JoiV2N\)Ji{V2N\) (8.26) 

The density of state p(A) is obtained by p(A) = V2NXp{X). The behavior 
is different from the complex case. The oscillations here are milder than the 
complex case. 
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9 Discussion 



We have explored the cross-over behavior near the origin for the density of 
state for hermitean and for real matrices made of blocks. Although the result 



has been known H, 20, 21, 22] for the one matrix model, our derivation, in 



particular our derivation using Kazakov's method, provides new and simple 
expressions. We have also extended the discussion to rings and lattices of 
matrices. We have proved the universality of this cross-over behavior to first 
order in the deviation from a Gaussian distributions, and conjecture that 
this universality should hold in general. 

We note also that in Kazakov's method we encounter expressions which 
are valid for a non-vanishing external source matrix A. Instead of letting all 
the eigenvalues of A go to zero, we may consider some specific choices of the 
eigenvalues of A. Let us, for instance, consider = cos{2TT'~f / N) . Then the 
contour integral UAit) in (|3.1CI|) depends upon these a^'s. However in the 
cross-over region, these a^'s may be neglected with respect to the variable u, 
which was scaled by a factor N in the integral. Thus such non-zero a-y's are 
irrelevant, and again we would obtain with them the same universal form in 
the large N limit. Unfortunately we have not been able to extend Kazakov's 
method to the case of a lattice of matrices. 

For the lattice of matrices, we could have considered also a representation 
of the Hamiltonian as a block matrix M of the form given in ( |2.1|) . Indeed, for 
a bi-partite lattice with an even number of sites, we could divide the lattice 
into two sub-lattices A and B such that nearest neighbors belong to different 
sub-lattices. The sub- lattice A contains sites labelled from 1 to N/2 and the 
sub-lattice B from N/2 -|- 1 to N. With such a labelling the Hamiltonian M for 
a particle hopping on this lattice would have a form like ( |2.1| ). However, the 
corresponding block matrix C would be sparse, containing many zero matrix 
elements. Indeed, the matrix C would be ^ ^ y; i-^- C would contain N'^/A 
elements. For a D-dimensional hypercubic lattice the number of bonds in 
the lattice is DN, where D is the spatial dimension. Then the ratio of the 
number of bonds to the number of matrix elements of C is given by 4D/N. 
When D = N/4, the model reduces to the one matrix model. The example 
A^ = 4 in the one-dimensional ring has been mentioned before. From this 
consideration, in the large D limit, if we keep D = N/4, we have a simple 
one matrix model of the form in ( |2.1| ) and the universality of the cross-over 
behavior near the center of the spectrum should hold. 
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As already remarked in one particular simple way to test for the 
universal oscillation studied here is to compare the height of the first peak 
in the density of states to the height of the second peak. According ( |7.17| ), 
the density of states is proportional to the universal function 

r{y) = y{Jl{y) + Jl{y)) (9.1) 

where y is the energy multiplied by some suitable constant. The positions of 
the peaks (and valleys) are determined by 

^ = Jl{y) - Jl{y) = (9.2) 
dy 

We find for example that the ratio of the height of the first peak to the height 
of the second peak is given by 1.218. (The ratio of the height of the first 
peak to the height of the first valley is 1.58.) 

A possible application is to the problem of a single particle propagating 
on a square lattice penetrated by random magnetic flux p3|, a problem that 



has recently attracted considerable attention Already the authors of 



(see figure 1 in this reference) noted that the density of states exhibits 
oscillations for finite A^. It is far from clear that our present work can be 
applied to this problem since, as we have just explained, the relevant matrix 
C in this random flux problem is sparse, with the ratio of non-zero matrix 
elements to the total number matrix elements given by 16/A^. Nevertheless, 
we note that a recent numerical study of the random flux problem by Avishai 



and Kohmoto found that the ratio of the height of the first peak to the 
height of the second peak was given by 1.25. We do not know whether the 
difference between 1.25 and 1.218 is real or due to numerical uncertainties. 
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APPENDIX A: INTEGRATION FOR GENERAL L 



We evaluate the following integral Dq for general L, which appeared in 



^0 ^=l 

6{t2+U + --- + tL-l) 



(27r)2 'L 



L 



(A.l) 



We can interpret the quadratic form in the square bracket as t'Ht with H 
the quantum Hamiltonian of a particle hopping on a ring (with a trivial 
constant site energy.) We diagonalize H and obtain for its eigenvalues the 
Bloch energies cos^j^. 

For example, in the case L = 4, we have 

JJ = ^ /^\4(Af-l)^-2Af f -2N{xl+xl+2xl) iki(xi+X4)+ik2{-xi+X4) 

° (27r)2^2^ J 
dkidk2dxidx2dxsdx4 

-2N {xl+xl+2xD 



.l)4(iV-l)g-27v| dx2dX3dXi5{2x^)e 



= ^(^)'^^-^^e-^ (A.2) 

where = ^(ti - t2 + h - U), X2 = ^(ti - ta), a^s = ^(t2 - ^4) and 
X4 = + ^2 + ^3 + ^4)- For general L, we obtain 

^0 = (f )^(^-^)e-#^(-^)^ n (l-cos(|fc))^i| (A.3) 
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Using the identity of ( p.29| ), we have 



Thus we get 



with 



n (l-cos^A:) = lnsin^(f^)2"-^ 

k=l,k^^ k=l 



L 



2 

2 



(A.4) 



Do = {j)'^^-%--^{^)'^D (A.5) 
D = ^ (A.6) 
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